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ABSTRACT 


Classical  concepts  of  Griffith-Irwin  fracture  mechanics  \ere  extended 
to  include  time  effects  and  hence  to  apply  to  crack  instability  under 
dynamic  loads.  Theoretical  considerations  of  the  stress  intensity 
histories  experienced  by  cracks  subjected  to  short  pulse  loads  suggested 
that  instability  requires  the  crack  tip  stress  intensity  to  exceed  the 
dynamic  toughness  for  a  minimum  time.  This  postulate  was  checked  by 
experiments.  Impact  techniques  were  used  to  produce  well-defined  stress 
pulses  in  specimens  of  an  epoxy,  4340  steel,  1018  steel,  and  6061-T6 
aluminum  containing  cracks  of  several  lengths  to  observe  crack  insta¬ 
bility  behavior.  The  experimental  results  were  not  well  described  by 
static  fracture  mechanics,  but  were  in  accord  with  the  minimum  time 
postulates  deduced  from  stress  intensity  histories.  The  understanding 
gained  from  the  research  program  promises  to  be  useful  in  assessing  the 
safety  of  structures  under  dynamic  loads,  in  characterizing  the  dynamic 
fracture  resistance  of  materials,  and  in  designing  equipment  and 
procedures  for  measuring  dynamic  fracture  toughness. 

The  relationship  between  dynamic  initiation  toughness  and  dynamic 
propagation  toughness  was  addressed.  The  temperature  histories  in 
materials  near  a  rapidly  propagating  crack  in  a  steel  plate  were 
measured  and  used  to  compute  the  heat  of  fracture  and  the  dynamic 
propagation  toughness.  Simultaneous  measurements  of  the  dynamic  stress 
intensity  were  made  using  the  shadow  optical  method  of  caustics  and 
high  speed  photography.  The  results  indicate  chat  more  energy  is 
delivered  to  the  running  crack  dp  chan  is  absorbed  there,  suggesting 
that  the  conventional  energy  balance  method  for  determining  dynamic 
propagation  toughness  may  lead  to  erroneous  results. 
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I  INTRODUCTION 


Structures  used  by  the  United  States  Air  Force  must  be  designed 
to  resist  catastrophic  fracture  when  subjected  to  dynamic  loads.  For 
example,  aircraft  components  may  experience  short  stress  pulses  from 
airborne  debris,  military  projectiles,  or  intense  bursts  of  laser  or 
x-ray  radiation.  Landing  gear  and  aircraft  retaining  cables  on  carrier 
ships  experience  dynamic  loads  at  the  end  of  each  flight.  To  ensure 
safe  design  of  such  structures,  a  knowledge  of  the  dynamic  fracture 
behavior  of  the  component  materials  is  necessary. 

This  final  technical  report  summarizes  the  results  of  a  four-year 
research  program  aimed  at  improving  our  understanding  of  dynamic 
fracture  behavior.  The  specific  research  objectives  are  given  in  the 
following  section.  Sections  III  and  IV  summarize  our  research  on 
dynamic  crack  instability  and  rapid  crack  propagation.  Section  V 
lists  the  publications  and  presentations  resulting  from  this  work,  and 
Section  VI  lists  the  project  personnel  and  describes  interactions  with 
other  scientists  nationally  and  internationally. 


4 


II  RESEARCH  OBJECTIVES 


This  research  program  had  two  objectives.  The  first  was  to  extend 
classical  Griff ith-Irwin  fracture  mechanics  by  accounting  for  time 
effects  to  obtain  instability  criteria  for  cracks  under  short  pulse 
loads.  The  second  objective  was  to  perform  critical  experiments  to 
measure  the  relation  between  energy  delivered  to  and  energy  absorbed 
at  the  tip  of  a  rapidly  running  crack.  Specific  research  tasks  were 
as  follows: 

(1)  Produce  arrays  of  well-defined  internal  macrocracks  of 
given  sizes,  shapes,  and  spacing  in  a  polymeric  model 
material . 

(2)  Load  the  cracks  with  well-defined  rectangular  stress  pulses 
and  thereby  generate  a  data  base  for  dynamic  crack 
instability. 

(3)  Use  the  data  to  evaluate  a  criterion  for  dynamic  crack 
instability  based  on  stress  intensity  considerations. 

(4)  Generate  dynamic  crack  instability  data  in  several  structural 
alloys  using  single-edge-nocehed  strip  specimens  and  varying 
stress  pulse  duration,  and  compare  results  with  the  dynamic 
crack  instability  criterion. 

IS)  Recommend  a  simplified  procedure  for  measuring  dynamic 

fracture  toughness  that  can  become  a  standard,  widely  used, 
laboratory  test. 

(6)  Begin  to  Investigate  the  relationship  between  dynamic  initia¬ 
tion  toughness  and  dynamic  propagation  toughness  Kj^  by 
measuring  the  heat  generated  by  a  rapidly  running  crack. 

In  successfully  completing  all  these  tasks,  we  have  made 
considerable  progress  in  understanding  the  dynamic  fracture  behavior 
of  materials.  Major  accompl i  aliments  include  the  following: 

•  formulation  of  stress  intensity  history  concepts  that  can 
explain  the  instability  behavior  of  cracks  loaded  by  short¬ 
lived  stress  pulses. 
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III  DYNAMIC  CRACK  INSTABILITY 

Background  and  Theoretical  Considerations 

According  to  classical  static  fracture  mechanics  (Paris  and  Sih, 

1964,  for  example),  the  stress  a  necessary  to  cause  crack  instability 

* 

varies  inversely  as  the  square  root  of  the  crack  length  a.  The  usual 
instability  criterion  compares  the  crack  tip  stress  intensity  K  with 
the  fracture  toughness  K^  such  chat  instability  occurs  when 

K  «  f(o,a1/2)  >  Kc 

The  fracture  toughness  is  a  general  measure  of  the  resistance  of  a 
material  to  unstable  crack  growth.  Under  plane  strain  conditions, 

K  (■  KIt>  is  geometry-independent  but  may  vary  with  strain  rate  and 
temperature. 

When  the  load  is  applied  very  rapidly,  the  material  may  exhibit  a 

different  fracture  resistance  K  ,♦  and  if  the  load  duration  is  short  and 

d 

comparable  to  the  time  required  for  a  stress  wave  to  run  the  length  of 
the  crack,  the  crack  tip  stress  intensity  may  vary  with  time.  This 
work  focused  on  the  time-dependence  of  the  stress  intensity  and  its 
effect  on  the  instability  condition.  As  for  the  material-related  part, 
a  constant  dynamic  fracture  toughne:  a  value  was  assumed  to  exist  and  to 
characterize  the  material  for  the  strain  race  considered. 

When  a  stress  wave  strikes  a  crack,  a  complicated  pattern  of 
diffracted  waves  is  generated  and  produces  initial  oscillations  in  the 
crack  tip  stress  field.  The  i lmo-dependent  crack  tip  stress  intensities 
produced  by  impinging  step  function  loads  have  been  computed  by  Sih, 

Kmbly.  and  Raver*  (1968,  1969,  1971,  1972)  and  Arhenbach  (1970,  1972)  and 
others.  These  results,  shown  in  Figure  I,  shot.*  that  the  dynamic  stress 
intensity  K^n  rises  sharply  with  time  (K^n  *  t^‘),  overshoots  the 

For  simplicity  of  discussion,  we  consider  here  only  Mode  1  (opening  mode) 
uniform  tensile  h  .filing  of  a  sharp  crack  in  a  perfectly  elastic  material 
and  neglect  edge  effects. 
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equivalent  static  stress  intensity  K  by  a  considerable  amount,  and 
then  reaches  a  constant  static  value  after  several  damped  oscillations. 

At  early  times,  is  independent  of  crack  length,  depending  only  on 

time  and  on  the  amplitude  of  the  stress  pulse.  This  is  because  the 
effective  crack  length  that  contributes  to  the  crack  tip  stress  intensity 
increases  linearly  with  time  as  the  stress  wave  propagates  along  its 
length  and  loads  it,  A  single-valued  parameter  that  characterizes  this 
complicated  stress  field  is  needed  to  allow  comparisons  with  the  material 
property  and  to  obtain  an  instability  criterion. 

Sih  (1968,  1973)  and  Achenbach  (1970,  1972)  speculated  that  dynam¬ 
ically  loaded  cracks  may  become  unstable  at  stress  levels  lower  than  those 
required  for  cracks  under  equivalent  static  loads  because  K^n  briefly 

S  C  cl  C 

exceeds  K  ‘  .  These  authors  thus  suggested  that  the  maximum  value  of 
K^n  determines  crack  instability. 

An  experiment  reported  by  Shoe key  and  Curran  (1973),  however,  suggested 

that  neither  the  static  fracture  mechanics  nor  the  K^n  criterion  could 

max 

explain  dynamic  crack  instability.  In  this  experiment  a  plate  specimen 
of  polycarbonate  containing  a  distribution  of  internal  penny-shaped 
cracks  was  impacted  by  a  flying  plate  to  produce  a  tensile  pulse  of 
known  amplitude  and  duration.  Cracks  above  a  certain  size  propagated, 
whereas  smaller  cracks  remained  stable,  and  the  dynamic  fracture  tough¬ 
ness  was  determined  from  the  critical  crack  size,  the  amplitude  of  the 
tensile  pulse,  and  the  static  fracture  mechanics  formula.  However,  when 
these  results  were  used  to  calculate  the  crack  instability  behavior  ob¬ 
served  in  a  similar  experiment  with  a  lower  amplitude  pulse,  a  discrepancy 

was  evident,  demonstrating  the  inadequacy  of  both  che  static  fracture 

d  vn 

mechanics  criterion  and  the  K  *  criterion  to  describe  track  instability 

max 

under  short  pulse  loads. 

Kalthoff  and  Shockey  (1977)  explained  these  results  by  analyzing  the 
early-time  stress  intensity  histories  experienced  by  cracks  of  different 
lengths  under  stress  pulses  of  different  durations  and  proposed  that, 
for  instability  to  occur,  must  exceed  K^,  the  dynamic  toughness, 

for  a  certain  minimum  time.  These  authors  postulated  that  static  fracture 
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mechanics  can  predict  instability  behavior  when  the  ratio  of  pulse 
duration  T0  to  crack  length  a  is  such  that  c^T0/a  >20,  where  is  the 
longitudinal  wave  speed.  Short  cracks  that  can  be  traversed  several 
times  by  a  wave  during  the  life  of  the  stress  pulse  experience  essenti¬ 
ally  constant  stress  intensity  histories  very  similar  to  those  produced 
by  static  loads;  hence  their  instability  response  is  governed  by  the 
classical  static  fracture  mechanics  formulas. 

Because  larger  cracks  experience  more  of  the  transient  aspects  of 
thj  stress  intensity  history,  the  dynamic  instability  curve  in  stress/ 
crack  length  space  lies  above  the  static  curve,  0  °Va.  For  cracks  above 
a  certain  length,  postulated  to  be  ciTQ/a<1.6,  constant  behavior  is 
predicted,  since  all  cracks  greater  than  this  size  experience  identical 
stress  intensity  histories. 

These  instability  deductions  for  penny-shaped  cracks  are  diagrammed 
in  Figure  2  as  a  surface  in  stress  amplitude/stress  duration/crack  length 
space.  Unstable  crack  growth  is  expected  for  conditions  of  0,  T0,  and 
a  that  lie  above  the  surface.  For  short  crack  lengths  and  long  stress 
durations,  the  instability  surface  is  defined  by  the  static  criterion. 
However,  for  long  crack  lengths  and  short  stress  durations,  time,  and  not 
crack  length,  is  important.  In  this  regime,  the  instability  stress  is  a 
function  only  of  the  stress  duration,  so  all  cracks  longer  than  a  certain 
length  exhibit  similar  behavior.  To  check  these  hypotheses,  we  applied 
stress  wave  loads  to  cracks  in  an  epoxy,  in  two  steels,  and  in  an 
aluminum  alloy  and  observed  their  instability  response.  The  results  are 
described  in  the  following  paragraphs. 

Experiments  and  Results 

Epoxy 

A  data  base  for  dynamic  crack  instability  was  generated  by  casting 

A 

arrays  of  appropriately  sized  Mylar  or  Kapton  disks  to  serve  as 


* 

Trademarks,  E.  I.  du  Pont  de  Nemours  &  Co.,  Wilmington,  Delaware  19898. 
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FIGURE  2  CRACK  INSTABILITY  SURFACE  ACCORDING  TO  DYNAMIC 
MINIMUM-TIME  STRESS  INTENSITY  CRITERION 
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internal  penny-shaped  cracks  in  epoxy  (Figure  3a)  and  subjecting  the 
specimens  to  rectangular,  2.04-ps-duration  stress  pulses  (Shockey  and 
Erlich,  in  preparation).  Six  plate  impact  specimens  of  Epon  Resin  815 
epoxy*  were  cast  from  a  single  batch  of  epoxy  in  an  effort  to  minimize 
specimen-to-specimen  variation  in  mechanical  properties.  The  specimens 
were  6.35  mm  thick  by  63.5  mm  in  diameter,  and  each  contained  six 
judiciously  spaced,  50-pm-thick  circular  Mylar  disks  having  radii  of 
6.35,  3.18,  1.59,  0.79,  and  0.40  (two  disks)  mm.  A  light  gas  gun  wa3 
used  to  accelerate  a  flyer  plate  against  the  specimen  plates  to  achieve 
rectangular  stress  pulses  of  2.04-ys  duration.  Impact  velocities  ranging 
from  0.0110  to  0.0275  mm/us,  corresponding  to  peak  tensile  stresses** 
ranging  from  14.9  to  37.3  MPa,  were  chosen  in  an  attempt  to  make  differ¬ 
ent  crack  sizes  grow  in  each  specimen.  (Examples  of  the  crack  growth 
induced  from  Kapton  disks  are  shown  in  Figure  3(b)  and  (c).) 

Figure  4  presents  the  data  in  stress-crack  radius  space  at  constant 
stress  duration  T0  =  2.04  ys.  The  solid  points,  indicating  crack  growth, 
and  the  open  circles,  denoting  no  crack  growth,  define  the  regimes  of 
crack  instability  and  stability,  respectively.  The  boundary  between  the 
two  regimes  represents  the  threshold  conditions  for  crack  instability  and 
defines  the  curve  required  to  evaluate  proposed  instability  criteria. 

Bev-ause  some  scatter  exists  in  the  data,  comparisons  of  proposed 
criteria  must  be  made  with  a  band  rather  than  a  sharp  curve. +  Neverthe¬ 
less,  the  band  provides  a  reliable  indication  of  the  shape  and  average 
location  of  the  instability  curve  and  thus  allows  discriminative  evaluations 
of  the  proposed  criteria. 

Epon  Resin  815,  Trademark,  Shell  Corporation,  Houston,  Texas. 

A  one-dimensional  elastic  w«f*  propagation  code,  SWAP,  computed  the 
durations  and  amplitudes  of  u.,-  peak  tensile  stresses  experienced  by 
the  individual  specimens. 

Scatter  is  attributable  in  part  to  the  relatively  large  inherent  varia¬ 
bility  in  mechanical  properties,  as  has  been  found  in  static  fracture 
toughness  data,  *  1.1  i  0.46  MPa 
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(a) 


(c) 


MP-G130-2 

FIGURE  3  CRACK  GROWTH  ABOUT  EMBEDDFD  KAPTON  DISKS 
RESULTING  FROM  APPLICATION  OF  A  SHDRT 
TENSILE  PULSF. 
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FIGURE  4  INSTABILITY  RESPONSE  OF  CRACKS  IN  EPOXY 

SUBJECTED  TO  2.04-/is-DURATION  STRESS  PULSES 
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Structural  Materials 

Cracks  in  single-edge-notched  strip  specimens  (1016  x  88.9  x  9.5  mm) 
of  4340  steel,  1018  steel,  and  6061-T651  aluminum  were  subjected  to  40- 
ps-long  haversine-shaped  tensile  pulses  of  various  amplitudes.  Pulse 
loads  were  applied  by  pneumatically  accelerating  a  cylindrical  projectile 
along  the  specimen  and  against  a  massive  block  attached  to  one  end  of  the 
specimen.  This  caused  a  tensile  wave,  whose  profile  was  measured  by  a 
strain  gage,  to  run  back  to  the  crack  location.  The  amplitudes  of 
successive  pulses  were  increased  in  small  steps  until  incremental  crack 
growth  could  be  detected  by  replicating  the  crack  tip  on  the  side  of 
the  specimen. 

The  results  are  shown  in  Figure  5.  For  all  three  materials,  the 
stress  necessary  for  crack  instability  decreased  strongly  with  crack 
length  (normalized  by  the  specimen  width  w  in  Figure  5)  at  short 
crack  lengths,  but  tended  to  reach  a  constant  value  for  crack  lengths 
greater  than  about  60%  of  the  specimen  width.  The  instability  curve 
predicted  by  the  static  criterion  (dashed  line)  is  shown  for  comparison 
and  is  dicussed  below.  Thus  similar  trends  in  the  variation  of  critical 
stress  with  eiaek  length  were  observed  in  epoxy,  4340  steel,  1018  steel, 
and  aluminum. 

Evaluation  of  Instability  Criteria 

Figure  6  compares  the  experimentally  observed  instability  curve  for 
epoxy  with  curves  predicted  by  several  proposed  criteria.  All  criteria 
predict  well  the  observed  behavior  at  small  crack  lengths  (a0  >0.3  mm, 
ciT0/a0  >20).  In  this  region,  the  stress  wave  has  time  to  run  the  length 
of  the  crack  several  times  before  the  pulse  vanishes,  so  the  crack  tip 
stress  field  approaches  equilibrium  and  static  conditions  apply.  For 
longer  crack  lengths  where  crack  tip  stress  conditions  are  further  from 
equilibrium,  the  predictions  of  all  criteria  are  progressively  worse. 

The  continuous  decrease  in  critical  stress  with  increasing  crack  length 
predicted  by  static  fracture  mechanics  disagrees  with  the  constant 
critical  stress  observed  for  cracks  larger  than  about  2  mm.  The  maximum 
dynamic  stress  intensity  criterion  predicts  constant  behavior  for 


STRESS  PULSE  AMPLITUDE, 


MA-61 30-240 


FIGURE  5  OBSERVED  INSTABILITY  BEHAVIOR  OF  CRACKS 
IN  THREE  STRUCTURAL  MATERIALS  LOADED 
BY  A  AO-m  STRESS  PULSE  AND  COMPARISON 
WITH  BEHAVIOR  PREDICTED  BY  STATIC 
FRACTURE  MECHANICS 
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CALCULATED  PEAK  TENSILE  STRESS,  an  —  MPa 


DIMENSIONLESS  CRACK  LENGTHS, 
CiVa0  =  5.3  mm/aQ 


13.3  6.71  3.33  1.67  0.835 


MA-6130-26A 


FIGURE  6  COMPARISON  OF  EXPERIMENTALLY  DETERMINED 
DYNAMIC  CRACK  INSTABILITY  CURVE 
WITH  CURVES  PREDICTED  BY  PROPOSED 
INSTABILITY  CRITERIA 
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tracks  greater  than  2.5  mm,. but  the  magnitude  is  significantly  less  than 
observed.  Moreover,  the  magnitude  of  the  critical  stress  predicted  for 
shorter  cracks  is  also  less  than  observed. 

The  data  are  best  described  by  the  minimum  time  criterion,  which 
accurately  predicts  the  shape  of  the  instability  curve,  but  slightly 
underpredicts  the  critical  stress  in  the  long  crack  length  region. 

More  consideration  is  needed  to  explain  the  observed  discrepancy  at 
long  crack  lengths. 

For  the  three  structural  materials,  results  similar  to  those  in 
epoxy  were  obtained.  Figure  5  shows  that  the  amplitudes  of  the  critical 
pulses  decreased  with  increasing  crack  length  at  a  rate  much  less  than 
predicted  b>  the  classical  fracture  mechanics  formula.  Furthermore, 
the  maximum  value  of  the  dynamic  stress  intensity  (as  calculated  by  a 
dynamic  finite  element  method  (Homma  and  Shockey,  in  preparation) 
varied  with  crack  length  for  the  two  steels,  and  thus  does  not 

appear  to  be  a  useful  dynamic  crack  instability  criterion. 

The  structural  alloy  results  are,  however,  in  qualitative  accord 
with  the  minimum  time  instability  criterion.  The  instability  curves 
at  constant  stress  duration  fall  rapidly  with  crack  length  at  short 
crack  lengths  and  then  appear  to  approach  a  constant  stress  value  at 
longer  crack  lengths  (compare  Figures  3  and  5).  Thus,  the  observed 
shape  of  the  instability  curve  agrees  with  the  predicted  shape,  and 
we  can  conclude  that  the  crack  instability  criterion  of  classical  static 
fracture  mechanics  can  be  modified  to  apply  to  high-rate,  short-duration 
load  situations. 

Subsequent  research  efforts  were  directed  at  obtaining  crack 
instability  data  at  both  shorter  and  longer  pulse  durations  and,  thereby, 
mapping  the  instability  surface  in  crack-length,  stress,  pulse  duration 
space.  The  dynamic  instability  behavior  of  cracks  in  4340  steel  strip 


*For  cjT0/n0  >  2.2,  the  stress  for  instability  stays  constant  although  } 

the  crack  length  increases,  because  the  stress  pulse  vanishes  before  the 
loading  wave  can  run  the  length  of  the  crack. 
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specimens  under  18-ys  and  8C-ys  stress  pulses  was  ascertained  by  the 
experimental  procedure  described  earlier.  The  results,  shown  in  Figure  7, 
support  the  minimum  time  criterion.  The  larger  pulse  durations  allow 
longer  cracks  to  follow  the  parabolic  static  instability  curve,  whereas 
the  short  18-ys  duration  pulses  produced  the  flat  dynamic  response  curve 
for  crack  lengths  as  short  as  0.2. 

These  results,  taken  together  with  the  earlier  results  for  40-ys 
pulses  (Figure  5a),  can  be  used  to  construct  the  instability  surface 
for  4340  steel,  as  shown  in  Figure  8.  Comparison  of  this  surface  with 
the  surface  predicted  from  theoretical  considerations  of  the  stress 
intensity  histories.  Figure  2,  shows  the  similarity  and  thus  further 
verifies  the  theoretical  concepts. 

Routine  Test  Procedure 

The  understanding  of  dynamic  crack  instability  gained  from  this 
research  program  can  be  used  to  design  a  testing  apparatus  and  outline 
a  testing  procedure  for  routine  evaluation  of  the  dynamic  fracture 
properties  of  materials.  To  be  widely  used,  the  testing  procedure  must 
be  simple,  and  the  test  facility,  specimens,  and  instrumentation  must 
be  inexpensive. 

We  suggest  modifying  a  commercially  available  dynamic  testing 
machine,  such  as  a  standard  Charpy  pendulum  or  dropweight  tester,  to 
produce  and  record  tensile  pulses  having  durations  ranging  from  10  to 
90  ys  and  amplitudes  up  to  500  MPa.  The  durations  require  a  series  of 
striker  plates  of  various  dimensions;  the  stress  amplitudes  will  require 
higher  Chan  normal  striking  velocities,  which  may  be  achieved  with 
springs.  The  production  of  clean  and  isolated  stress  pulses  requires 
simultaneous  impact  of  the  striker  plates.  When  the  striker  plate  is 
short,  the  shape  of  the  pulse  depends  strongly  on  the  lateral  striker 
plate  dimensions.  The  design  of  Che  striker  assemblies,  guided  by 
intuition  in  our  previous  work,  can  be  accomplished  by  two-dimensional 
finite  difference  wave  propagation  codes.  The  abilities  of  the  strikers 
to  produce  the  desired  stress  pulses  can  be  tested  by  mounting  the 
strikers  on  the  impact  machine  and  using  strategically  located  strain 
gages  to  record  the  stress  pulses  resulting  from  the  impact. 
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FIGURE  7  EFFECT  OF  PULSE  DURATION 
ON  INSTABILITY  STRESS  FOR 
CRACKS  IN  4340  STEEL 
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Amplitudes  of  subsequent  tensile  stress  pulses  must  be  less  than 
the  initial  one.  In  the  present  work,  adequate  damping  of  subsequent 
tensile  pulses  was  achieved  by  gradually  decelerating  striker  plates 
and  projectiles  with  energy-absorbing  aluminum  honeycomb,  but  other 
reusable  and  tunable  energy-absorbing  systems  could  be  incorporated. 

Finally,  instrumentation,  if  needed,  may  be  attached  to  the  striker 
of  the  machine  to  avoid  having  to  instrument  each  specimen,  and  the 
striker  should  be  fixed  and  impacted  by  the  specimen.  These  ideas  will 
be  pursued  in  a  follow-on  contract  with  AFOSR. 


IV  DYNAMIC  CRACK  PROPAGATION 


Dynamic  fracture  toughness  can  mean  either  the  resistance  exhibited 
by  the  material  to  propagation  of  an  existing  crack  subjected  to  a  dynamic 
load  (i.e.,  a  stress  pulse)  or  the  resistance  exhibited  by  the  material 
to  a  rapidly  running  crack  under,  say,  a  static  load.  In  the  research 
described  in  the  preceding  section,  the  former  case,  dynamic  crack 
initiation,  was  investigated.  In  this  section  we  discuss  our  investigation 
of  the  latter  situation,  deducing  the  dynamic  fracture  toughness  associated 
with  a  rapidly  running  crack,  Kjq,  from  measurements  of  temperature  histories 
near  the  crack  path.  The  results  were  compared  with  the  energy  delivered 
to  the  crack  tip  as  measured  simultaneously  by  the  method  of  caustics. 

Background 

The  energy  absorbed  by  a  propagating  crack  is  absorbed  by  plastic 
flow  and  fracture  within  a  relatively  small  (compared  with  specimen 
dimensions)  volume  of  material  in  the  vicinity  of  the  crack  tip.  Calcu¬ 
lations  show  that,  in  metals,  the  energy  absorbed  in  creation  of  new 
surface  area  is  negligible  compared  with  the  plastic  energy;  therefore, 
the  fracture  toughness  associated  with  a  running  crack  is  determined  by 
the  plastic  energy  expended  at  the  crack  tip. 

The  plastic  energy  is  dissipated  as  heat,  which  flows  awav  from 
the  crack  plane  into  both  halves  of  the  specimen,  causing  small,  transient 
temperature  elevations.  If  the  temperature  history  of  material  near  a 
running  crack  can  be  measured,  the  heat  energy  and  hence  the  fracture 
toughness  of  the  running  crack  can  be  calculated  by  olving  tbe  appro¬ 
priate  heat  flow  equation.  Thus  the  results  will  yield  direct  measure¬ 
ments  of  dynamic  propagation  toughness  Kjd,  which  can  be  compared  with 
the  dynamic  initiation  toughness  measured  in  the  Impact  experiments 
in  this  program 

A  thermocouple  technique  was  used  to  measure  temperatures  associated 
with  rapidly  propagating  cracks  ard  to  obtain  fracture  energies  from  the 
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temperature  histories.  A  NATO  travel  giant  was  awarded  to  Dv.  D.  A.  Shockey, 
allowing  him  to  make  these  measurements  at  the  Institut  fur  Werkstoffmechanik 
in  Freiburg,  Germany,  where  he  was  able  to  make  simultaneous  measurements 
of  the  dynamic  stress  intensity  at  the  rapidly  moving  crack  tip  using 
unique  German  equipment  to  apply  the  method  of  caustics.  Expert  technical 
assistance  was  provided  by  J.  F.  Kalthoff,  W.  Klemm,  S.  Winkler,  and  A.  Wolf. 
This  effort,  which  was  highly  successful,  has  been  written  up  in  preparation 
for  submission  to  a  technical  journal  and  is  described  briefly  here. 

Specimens 

Six  compact  fracture  specimens  250  by  250  by  20  mm  were  machined 
from  a  high-strength  maraging  steel*  and  then  hardened  to  Rc  56  (yield 
strength  =  2100  MN  m  2)  by  heating  to  480°C  for  4  hours.  The  static 
fracture  toughness  Klc  was  measured  to  be  50  MN  m™3/2.  A  208-mm-long 
slot  having  a  45°  chevron  tip  ser\ad  as  the  starter  notch.  One  side  of 
each  specimen  was  lapped  and  polished  to  a  mirror  finish  in  reparation 
for  measuring  Kjyn  by  the  caustics  method.  Four  to  ten  30-pm-diameter 
constantan  wires  were  spot -welded  to  the  specimen  near  the  expected  crack 
path  to  monitor  temperature  histories  and  allow  subsequent  calculations 
of  Kjj).  Figure  9  shows  a  prepared  specimen  before  and  after  it  was  placed 
in  the  testing  machine.  The  specimens  were  loaded  in  a  hydraulically 
activated  mechanical  testing  machine  by  the  transverse  split  wedge  tech¬ 
nique.  A  range  of  crack  velocities  vas  achieved  by  causing  crack 
initiation  to  occur  under  different  specified  stress  intensities 

Dynamic  Stress  Intensity  Measurement 

The  energy  delivered  to  the  crack  tip  was  measured  by  the  shadow 
optical  method  of  caustics  (Kalthoff  et  al.,  1977).  In  this  method, 
light  incident  on  the  specimen  surface  is  reflected  into  a  camera.  A 
dark  shadow  spot  appears  at  the  crack  tip  because  the  surface  there  is 


*The  German  designation  for  the  steel  is  HFX  760.  It  has  a  nominal 
composition  of  18%  Ni,  9%  Co,  4.8%  Mo,  0.9%  Ti*  and  <0.03%  C  and 
corresponds  to  the  American  18  Ni  maraging  grade  300, 


deformed  inward  by  the  high  stresses  that  exist  at  the  crack  tip.  The 
size  of  the  shadow  spot  is  a  measure  of  the  crack  tip  stress  intensity 
K^yn  and  hence  a  measure  of  the  energy  available  at  the  crack  tip: 

Kjyn  =  H  F(v)  D5/2  (1) 

where  M  is  a  constant  that  depends  on  material  properties  and  the  experi¬ 
mental  arrangement,  F(v)  is  a  dimensionless  factor  that  corrects  for  the 
effect  of  crack  speed,  and  D  is  the  diameter  of  the  shadow  spot  normal 
to  the  fracture  plane. 

The  energy  delivered  to  the  crack  tip  as  the  crack  propagated  was 
determined  at  24  locations  by  photographing  the  shadow  spot  at  the 
running  crack  tip  every  few  microseconds  with  a  spark  camera. 

Heat  of  Fracture  Measurement 

The  energy  absorbed  at  the  tip  of  a  running  crack  was  determined 
from  the  temperature  histories  of  material  near  the  crack  plane  as 
measured  by  sensitive  thermocouples.  Four  to  ten  constantan  wires  30  pm 
in  diameter  were  spot-welded  to  Che  specimen  surfaces  at  locations  along 
the  expected  crack  path  and  slightly  (0.7  to  2.0  mm)  to  the  side  of  the 
expected  crack  path.  An  equal  number  of  constantan  wires  were  spot-welded 
to  a  plate  of  the  same  steel  to  serve  as  reference  thermocouples. 

The  eime-temperature-distance  data  were  used  to  solve  the  heat-flow 
equation  and  to  calculate  Che  heat  liberated  by  the  crack  at  various 
locations  along  its  path.  The  liberated  heat  represents  nearly  all  the 
energy  absorbed  by  the  crack  and  is  thus  a  measure  of  the  dynamic  fracture 
toughness  K^). 

To  obtain  approximate  values  of  the  heat  generated  by  the  propagating 
crack  Q  and  the  dynamic  fracture  toughness  Kjp,  we  used  the  measured 
maximum  temperature  values  Tmax  to  solve  the  one-dimensional  heat  flow 
equation:  2 

Q  a  £rn,aX  y  five  Pc  -  (1  -  u2)  ~~D-  (2) 

where  y  is  the  distance  between  the  fracture  surface  and  the  thermocouple, 
e  is  the  constant  2.718,  and  the  other  symbols  are  properties  of  HFX  760 
steel  given  below: 
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Density,  p 

8040  kg  m~3 

Specific  heat,  c 

477  J  kg~loC~1 

Poisson  ratio,  u 

0.32 

Elastic  modulus,  E 

183  x  1011  N  m' 

Results  and  Discussion 

Seven  experiments  were  performed  at  crack  velocities  ranging  from 
500  to  1200  m  s-^.  A  typical  series  of  shadow  optical  photographs  of 
a  propagating  crack  tip  is  shown  in  Figure  10.  The  diameter  of  the 
shadow  spot  in  each  frame  is  carefully  measured,  and  the  dynamic  stress 
intensity  K^n  at  the  crack  tip  at  that  instant  of  time  was  calculated 
from  equation  (1). 

A  record  of  the  temperature  history  obtained  by  a  thermocouple  is 
given  in  Figure  11.  No  change  in  temperature  is  detected  initially  after 
triggering  because  the  heat  generated  by  the  crack  requires  time  to  flow 
to  the  thermocouple  location.  When  the  heat  wave  arrives,  the  temperature 
rises  sharply  to  a  maximum  (typically  several  tenths  of  a  centigrade 
degree)  and  then  falls  off  slowly  as  the  heat  pulse  continues  to  dissipate 
in  the  specimen. 

The  first  five  experiments  produced  non-steady  crack  velocities  and 
substantial  scatter  in  the  Kjp  and  data.  We  suspected  that  residual 

stresses  were  responsible,*  so  the  final  two  experiments  (17  and  18) 
were  performed  on  specimen  plates  that  were  given  a  stress  relief 
anneal.  Spark  camera  photographs  of  the  propagating  crack  in  the  stress- 
relieved  specimens  showed  two  steady-state  crack  velocities,  with  the 
initial  velocity  being  slightly  greater  than  the  second,  as  is  typical 
of  compact  specimens.  Furthermore,  the  scatter  in  Kj^n  values  was 
markedly  reduced.  These  results  show  the  important  influence  of  reaidual 
stresses  on  dynamic  toughness  evaluation. 

*The  heat  treatment  and  subsequent  quench  given  to  the  specimens  to  obtain 
the  desired  mechanical  properties  caused  the  plates  to  warp  slightly,  which 
suggested  the  existence  of  residual  stresses. 
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Figure  12  shows  the  variation  of  Kj^n  and  Kjp  with  crack  length  for 
the  two  stress-relieved  specimens.  The  Kjq  values  are  significantly 
smaller  (10%  to  50%)  than  the  values,  indicating  that  considerably 
more  energy  is  delivered  to  the  crack  tip  than  is  absorbed  there.  How¬ 
ever,  these  results  must  be  regarded  as  preliminary  because  of  uncertain¬ 
ties  in  the  residual  stress  state  of  the  specimens  and  because  an 

k 

approximate  heat-flow  analysis  was  used.  Additional  work  is  necessary 
to  more  accurately  assess  the  relationship  between  the  energy  delivered 
to  and  the  energy  absorbed  by  rapidly  propagating  cracks. 


£ 

The  temperature  histories  calculated  assuming  one-dimensional  heat 
flow  did  not  replicate  the  measured  histories,  and  so  the  Kjp  values 
obtained  from  equation  (2)  are  in  some  error.  A  three-dimensional 
analysis  of  heat  flow  in  the  specimen  geometry  is  needed  to  obtain 
more  accurate  values  of  K^. 
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FIGURE  12  COMPARISON  OF  ENERGY  DELIVERED  TO  AND  ENERGY 
ABSORBED  BY  A  RAPIDLY  PROPAGATING  CRACK 
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"Criterion  for  Crack  Instability  Under  Short  Pulse  Loads," 

Advances  in  Fracture  Research,  D.  Francois  et  al.,  Eds.  (Oxford 
and  Pergamon  Press,  New  York)  pp.  415-423. 

H.  Hotnma  and  D.  A.  Shockey,  "Response  of  Cracks  in  Structural 
Materials  to  Short  Pulse  Loads — Experimental  Results,"  to  be 
submitted  to  the  Journal  of  the  Mechanics  and  Physics  of  Solids. 

K.  Homma  and  D.  A.  Shockey,  "Response  of  Cracks  in  Structural 
Materials  to  Short  Pulse  Loads — Calculations  and  Analyses,"  to 

be  submitted  to  the  Journal  of  the  Mechanics  and  Physics  of  Solids. 

D.  A,  Shockey,  and  D.  C.  Erlich,  "Instability  Behavior  of  Cracks 
Under  Short-Duration  Loads,"  to  be  submitted  to  the  International 
Journal  of  Fracture  Mechanics. 

D.  A.  Shockey,  J.  F.  Kalthoff,  W.  Klemm,  and  S.  Winkler,  "Simultane¬ 
ous  Measurements  of  Stress  Intensity  and  Toughness  for  Fast  Running 
Cracks  in  Steel,"  to  be  submitted  to  Experimental  Mechanic^. 


Presentations 

D.  C.  Erlich  and  D.  A.  Shockey,  "Instability  Conditions  for  Cracks 
Under  Short-Duration  Pulse  Loads,"  Topical  Conference  of  Shock 
Waves  in  Condensed  Matter,  meeting  of  the  American  Physical  Society, 
Washington  State  University,  Pullman,  Washington,  June  11-13,  1979. 

D.  A.  Shockey,  "Instability  Conditions  for  Cracks  Loaded  by  Short 
Stress  Pulses,"  Poulter  Seminar,  SRI  International,  Menlo  Park,  CA, 
December  12,  1979. 

D.  A.  Shockey,  "Dynamic.  Crack  Instability,"  Institut  CERAC,  Ecublens, 
Switzerland,  May  19,  1980. 


D,  A.  Shockey,  "Dynamic  Crack  Instability,"  Institut  fur  Werkstoff- 
mechanik,  Freiburg,  Germany,  May  21,  1980. 

D.  A.  Shockey,  "Simultaneous  Measurements  of  Stress  Intensity  and 
Toughness  for  Fast  Running  Cracks  in  Steel,"  Poulter  Laboratory 
Seminar,  SRI  International,  Menlo  Park,  CA,  June  1980. 

D.  A.  Shockey,  "Criterion  for  Crack  Instability  Under  Short  Pulse 
Loads,"  Symposium  on  Crack  Formation  and  Propagation,  sponsored 
by  Polish  Academy  of  Science,  International  Union  of  Theoretical 
and  Applied  Mechanics  and  the  International  Center  for  Mechanical 
Sciences,  Tuczno,  Poland,  March  23-27,  1981. 

D.  A.  Shockey,  "Criterion  for  Crack  Instability  Under  Short  Pulse 
Loads,"  Fifth  International  Conference  on  Fracture  (ICF5)  Cannes, 
France,  March  29-April  3,  1981. 


VI  PERSONNEL  AND  ACTIVITIES 


The  following  professional  personnel  were  associated  with  this 
research  effort. 

D.  A.  Shockey,  principal  investigator 

D.  C.  Erlich 

H .  Homma 

J.  F.  Kalthoff 

W,  Kletnm 

S.  Winkler 

Five  publications  in  scientific  journals  will  result  from  this 
program.  Air  Force  personnel  have  been  kept  abreast  of  technical 
developments  through  the  annual  reports,  and  contact  was  maintained 
throughout  the  program  with  Dr.  A.  H.  Rosenstein  of  the  Air  Force 
Office  of  Scientific  Research.  Dr.  Theodore  Nicholas,  Metals  Behavior 
Branch  of  the  Air  Force  Materials  Laboratory  at  Wright-Patterson  Air 
Force  Base,  was  periodically  consulted  by  telephone,  and  he  visited 
Dr,  Shockey  at  SRI  International  in  November  1980  to  obtain  a  complete, 
first-hand  account  of  the  accomplishments  in  the  program. 

This  program  has  attracted  international  attention.  The  principal 
investigator,  D.  A.  Shockey  has  been  invited  to  a  special  dynamic 
fracture  conference  in  Poland  in  1981  to  present  and  discuss  the  results. 

A  paper  describing  the  results  of  die  dynamic  crack  instability 
experiments  has  been  accepted  for  presentation  and  publication  at  the 
Fifth  International  Conference  on  Fracture  to  be  held  March  1981  in 
Cannes,  France. 

Furthermore,  the  North  Atlantic  Treaty  Organization  granted  travel 
funds  and  living  expenses  to  the  principal  investigator  to  allow  him  to 
use  a  unique  optical  arrangement  in  Germany  in  his  dynamic  fracture 
experiments.  These  experiments  were  highly  successful  and  are  reported 
in  Section  IV  of  this  report.  A  paper  describing  this  work  and  the  results 
has  been  prepared  and  will  be  submitted  for  publication  in  a  technical 
journal. 


34 


REFERENCES 


Achenbach,  J.  D.  (1970).  Brittle  and  ductile  extension  of  a  finite 

crack  by  a  horizontally  polarized  shear  wave.  Int.  J.  Eng.  Sci., 

8,  947-966.  '  '  . . . 

Achenbach,  J.  D.  (1972).  Dynamic  effects  in  brittle  fracture.  In 
Mechanics  Today,  S.  Nemat-Nasser  (Ed.),  Vol.  1,  Pergamon. 

Embley,  G.  T.,  and  G.  C.  Sih  (1971).  Proc.  12th  Midwestern  Mechanics 
Conference,  6,  473,  See  G.  C.  Sih  (1973).  Handbook  of  Stress 
Intensity  Factors,  Institute  of  Fracture  and  Solid  Mechanics, 

Lehigh  University,  Bethlehem,  PA.  p.  8, 5. 2-1. 

Homma,  H.  and  D.  A.  Shockey,  Response  of  cracks  in  structural  materials 
to  short  pulse  loads,  to  be  submitted  to  the  J.  Mech.  Phys.  Solids. 

Kalthoff,  J.  F.,  J.  Beinert  and  S.  Winkler,  (1977).  "Measurements  of 
Dynamic  Stress  Intensity  Factors  for  Fast  Running  and  Arresting 
Cracks  in  Double-Cantilever-Beam  Specimens,"  in  ASTM  STP  627,  Fast 
Fracture  and  Crack  Arrest,  The  American  Society  for  Testing  and 
Materials,  Philadelphia,  PA,  p.  161-176, 

Kalthoff,  J.  F.  and  D.  A.  Shockey  (1977).  Instability  of  cracks  under 
impulse  loads.  J.  Appl,  Phys.,  48  (3),  984-993. 

Paris,  P.  C.  and  G.  C.  Sih  (1964).  Stress  analysis  of  cracks.  In 
Fracture  Toughness  Testing  and  Its  Applications,  ASTM  Special 
Technical  Publication  No.  381. 

Ravera,  R.  J.  and  G.  C.  Sih  (1969).  Transient  analysis  of  stress  waves 
around  cracks  under  anti plane  strain.  J.  Acous.  Soc.  Amer. ,  47, 

875-880. 

Shockey,  D.  A.  and  D.  R.  Curran  (1973).  A  Method  for  Measuring  Kic  at 
Very  High  Strain  Rotes,  ASTM  Special  Technical  Publication  536. 

Fracture  Toughness  Testing,  the  ASTM,  Philadelphia,  PA,  p.  297-311. 

Shockey,  D.  A.  and  D.  C.  Erlich,  Instability  behavior  of  cracks  under 

short-duration  loads,  to  be  submitted  to  Internat.  Fracture  Mechanics 

Sih,  G.  C.  (3.968)  Some  elastodynamic  problems  of  cracks.  Int.  J,  Frac. 
Mech.,  4,  51-68. 

Sih,  G.  C.  (1973).  Handbook  t>f  Stress-Intensity  Factors,  Institute  of 
Fracture  and  Solid  Mechanics,  Lehigh  University,  Bethlehem,  PA. 

Sih,  G.  C. ,  G.  T.  Embley,  and  R.  J.  Ravera  (1972).  Impact  response  of 
a  plane  crack  in  extension.  Int.  J.  Solids  Structures,  8,  977-993. 


35 


DISTRIBUTION  LIST 


No.  Copies 

1.  AFOSR/NE  '  "  ‘ 

ATTN:  Dr.  A.  H.  Rosenstein  16 

Bldg.  #410 

Bolling  AFB 

Washington,  D.C.  20332 

2.  AFWAL/MLLM 

Attn:  Branch  Chief  1 

Wright-Pat terson  AFB,  Ohio  45433 

3.  AFWAL/MLLS 

Attn:  Branch  Chief  1 

Wright-Patterson  AFB,  Ohio  45433 

4.  AFWAL/MLLN 

Attn:  Ted  Nicholas  1 

Wright-Patterson  AFB,  Ohio  45433 

5.  AFWL 

Attn:  Joe  Renick  1 

Kir eland  AFB,  N.H.  87117 

6.  AFATL/DLJW 

Attn:  Maj.  Guy  Spitale  1 

Eglin  AFB,  Florida  32542 


36 


